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Modular inflation and the orthogonal axion as curvaton 
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We study a particular supersymmetric realization of the Peccei-Quinn symmetry which provides a 
suitable candidate for the curvaton field. The class of models considered also solves the /i problem, 
while generating the Peccei-Quinn scale dynamically. The curvaton candidate is a pseudo Nambu- 
Goldstone boson corresponding to an angular degree of freedom orthogonal to the axion field. Its 
order parameter increases substantially following a phase transition during inflation. This results 
in a drastic amplification of the curvaton perturbations. Consequently, the mechanism is able to 
accommodate low-scale inflation with Hubble parameter at the TeV scale. Hence, we investigate 
modular inflation using a string axion field as the inflaton with inflation scale determined by gravity 
mediated soft supersymmetry breaking. We find that modular inflation with the orthogonal axion 
as curvaton can indeed account for the observations for natural values of the parameters. 

PACS numbers; 98.80.Cq 



I. INTRODUCTION 

The last decade experienced a plethora of precise cos- 
mological observations. These observations have con- 
firmed the basic predictions of inflation rendering the 
latter an essential extension of the hot big bang model. 
Precision cosmology upgraded inflation model-building 
at the next level, that of the design of complex realistic 
models (in contrast to early single-field fine-tuned infla- 
tionary models). These models utilize and reflect the 
rich content of particle theory. A first such attempt was 
the well known hybrid inflation model [f-3], which cou- 
ples the inflaton field with the Higgs field of a grand 
unified theory (GUT). In this way, hybrid inflation dis- 
penses with a number of tuning problems, which plagued 
most simplistic single-field inflation models. 

In the same spirit, curvaton models [4-6] (see also [7]) 
use a second field to generate curvature perturbations. 
This curvaton field is not an ad hoc degree of freedom 
introduced by hand, but, instead, it is [8-10] a realistic 
field already present in simple extensions of the standard 
model (SM). In the context of the curvaton, inflation 
model-building is substantially liberated [II, 12] allowing 
for more realistic and less fine-tuned models. 

A particular advantage of inflation model-building in 
the context of the curvaton mechanism is that it is possi- 
ble to construct [II, 13] models with inflationary energy 
scale much lower than the GUT scale. In this spirit, we 
investigate, in this paper, the possibility of using a string 
axion field as the inflaton. 

The nature and origin of the inflaton field are still open 
questions in inflation model-building. Typically, what is 
required is a light field with suppressed interactions with 
other fields including those of the SM. By "light" , we 
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mean a field whose effective mass is smaller than the 
Hubble parameter i7» at the time when the cosmological 
scales exit the horizon during inflation. This guarantees 
that inflation lasts long enough to encompass the cosmo- 
logical scales. The interactions of this field have to be 
suppressed in order not to lift the flatness of the scalar 
potential along the inflationary trajectory. 

However, the case of slow-roll inflation (i.e. the case 
with inflaton mass <C -ff*) suffers from the fact that, typ- 
ically, supergravity (SUGRA) introduces [2, 14] correc- 
tions to the inflaton mass of order H^, during the infla- 
tionary period. To keep the inflaton mass under con- 
trol, one may use as inflaton a pseudo Nambu-Goldstone 
boson (PNGB) field, since the flatness of the potential 
of such a field is protected by a global U(I) symme- 
try. Promising such candidates are [15] the string ax- 
ions, which are the imaginary parts of string moduli fields 
with the flatness of their potential lifted only by (soft) 
supersymmetry (SUSY) breaking. This results in infla- 
ton masses of order H^. Hence, such modular inflation 
is of the fast-roll type [16]. Fast-roll inflation lasts only 
a limited number of e-foldings, which, however, can be 
enough to solve the horizon and flatness problems. 

The inflationary energy scale, in this model, is much 
lower than the GUT scale. As a result, the perturbations 
of the inflaton field are not sufficiently large to account 
for the observations. Consequently, a curvaton field is 
necessary to provide the observed curvature perturba- 
tion. However, even the curvaton cannot [17] generically 
help us to reduce the inflationary scale to energies low 
enough for modular inflation. This is possible only in 
certain curvaton models which amplify [13, 18] addition- 
ally the curvaton perturbations. We describe in detail 
such a model belonging to a class of SUSY realizations of 
the Peccei-Quinn (PQ) symmetry [19], which also solves 
the strong CP and fj, problems. We use as curvaton an 
angular degree of freedom orthogonal to the QCD axion. 

Our paper is structured as follows. In Sec. II, we 
present a brief outline of modular inflation. In Sec. HI, 
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wc analyze the amphfication mechanism for the curvature 
perturbations due to a PNGB curvaton with varying or- 
der parameter. In Sec. IV, we investigate whether it is 
possible to employ in the role of such a PNGB curva- 
ton an angular degree of freedom orthogonal to the QCD 
axion ("orthogonal axion") in SUSY theories. We show 
that the simplest constructions using two PQ superfields 
cannot work because the orthogonal axion is not appro- 
priately light during inflation. In Sec. V, we construct 
in detail an appropriate class of PQ models involving 
three SM singlet superfields, only two of which carry PQ 
charges. In Sec. VI, we study in detail the characteris- 
tics of the scalar potential in the above class of curvaton 
models. In Sec. VII, we focus on curvaton physics and de- 
rive a number of important constraints necessary for the 
model to be a successful curvaton model. In Sec. VIII, 
we quantify our findings in a concrete example of this 
class of models. We find that this concrete model can 
indeed work for natural values of the model parameters. 
Finally, in Sec. IX, wc discuss our results and present our 
conclusions. Throughout the paper, we use natural units, 
where c = h = 1 and the Newton's gravitational constant 
is G = 87rAfp ^ with Mp ~ 2.44 x 10^^ GeV being the re- 
duced Planck mass. 

II. MODULAR INFLATION 

String theory, in general, contains a number of mod- 
uli fields whose tree- level Kahler potential, in 4- 
dimensional effective SUGRA, is of the form 

/r=-M|5^1n[($i + $*)/Mp]. (1) 

i 

Hence, the Kahler potential is flat in the directions of the 
imaginary parts of the moduli Im($i). The same is true 
for the F-term scalar potential Vp, despite the fact that 
the superpotential may receive non-perturbative contri- 
butions (e.g. from gaugino condensation) of the form 

AI^($,) = A3exp(-A$,/Mp), (2) 

since Vp turns out to be independent of the phase of 
AW. The mass parameter A in Eq. (2) is usually taken 
to be the string scale and the Pi's are model-dependent 
coefficients of order unity. The Im(<I>i) fields are periodic 
(by modular invariance) 

Im($i) = lm($0 + 2wfi, (3) 

where fi ^ Mp. This is why the Im($j) fields are also 

called string axions. 

In compactified hetcrotic string theory, these axions 
correspond [20], in fact, to the massless modes of the 
second-rank antisymmetric tensor field B: 

B = b^^dx^ A dx" + biui^, dy'^ A dy^* , (4) 

where are the coordinates in the usual 4-dimensional 
space, y" the complex coordinates in the compactified 



extra-dimensional space and Lo^^^fj, harmonic (1,1) forms 
parametrizing the geometry of the compactified space 
[21]. The usual 4-dimcnsional components hp,i, of B cor- 
respond to the so-called model-independent axion, while 
the extra-dimensional components hj of B correspond to 
the so-called model-dependent axions. 

The flatness of the string axion potential is lifted only 
by soft SUSY breaking, which tilts the vacuum manifold 
by an amount determined by the SUSY-breaking scale. 
One can consider [15] that it is this potential that pro- 
vides the vacuum energy density necessary for inflation. 
Assuming gravity mediated soft SUSY breaking, the in- 
flationary potential V<t at the time when the cosmological 
scales exit the horizon is of intermediate scale: 

~ Y^ms/sMp ~ 10^°-^ GeV (5) 

for which iJ* ~ 7713/2, where 7713/2 ~ 1 TeV stands for the 
gravitino mass. The inflationary potential is of the form 

V{s) = V^-\my + ---, (6) 

where the ellipsis denotes terms which are expected to 

stabilize the potential at s ^ Mp with s being the canon- 
ically normalized string axion. Therefore, in the above 
formula, we have 

Kn ~ (to3/2-^p)^ and ~ m3/2. (7) 
This inflation model results in fast-roll inflation, where 

s = Si exp{FsAN) with Fs = ^ (^^/TT^c/9 - l) , 

Here, AA^ is the number of the elapsed e-foldings and Si 
the initial value of the inflaton field s. From the above, 
one can obtain the inflation potential A'' e-foldings before 
the end of inflation as 

ViN) :^V^{1- e-'^'^) . (9) 

Even though fast-roll, modular inflation keeps the 
Hubble parameter H rather rigid. Indeed, it can be easily 
shown that 

because c ~ 1 and s <i; Mp during inflation with e being 
one of the so-called slow-roll parameters defined as 



where the dot denotes derivative with respect to the cos- 
mic time. 

For modular inflation, the initial conditions for the in- 
flaton field are determined by the quantum fiuctuations 
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which send the field oH' the top of the potential hill. 
Hence, we expect that the initial value for the infiaton is 



"27 



(12) 



where ~ -/V^/v^Mp. 

Using the above and considering that the final value 
of s is close to its vacuum expectation value (VEV) 
svEV ~ Mp, we can estimate, through the use of Eq. (8), 
the total number of e- foldings as 



Ntnt =i ^ In 



Fs \m3/2 
where we took into account that 



Hjr 



m3/2- 



(13) 



(14) 



Inflation at such low energy scale as in Eq. (7) can 
provide the required amplitude for the curvature pertur- 
bations only through the use of a special kind of curvaton 
field, whose perturbations are amplified during inflation. 
In the following, we describe the mechanism for such am- 
plification. 



III. AMPLIFYING THE CURVATON 
PERTURBATIONS 

We discuss here the case of an axion-like curvaton, i.e. 
a PNGB. Examples of such a curvaton can be found in 
Refs. [9, 10]. However, in contrast to those works, we 
consider a PNGB curvaton whose order parameter has 
[13; 18] a different (larger) expectation value in the vac- 
uum than during inflation and, in particular, when the 
cosmological scales exit the horizon. Thus, the potential 
for the real canonically normalized curvaton field a is 



1 — cos 

1 



V{\a\«v)::.^mW 



(15) 



where v = v{t) is the order parameter of the PNGB (de- 
termined by the values of the radial fields in the model) 
with t being the cosmic time and rha- = rha-iv) is the mass 
of the curvaton at a given moment. In the true vacuum, 
we have v = vo and = with vo and m„ being the 
order parameter and the mass of the PNGB curvaton in 
the vacuum respectively. 

Note that, in principle, the PNGB does not need to 
have an exact sinusoidal potential. Instead, one could 
substitute [1 — cos((t/w)] by /(cr/w), where / is a peri- 
odic function of period 27r with a global minimum at the 
origin and /(O) = 0. Then, generically, the second line 
in Eq. (15) continues to be valid if a is close enough to 
the global minimum [23]. 



A. The amplification factor 

In this section, we will demonstrate that the curvaton 
perturbations can be amplified by the non-trivial evolu- 
tion of its order parameter v if the curvaton is a PNGB. 
This mechanism was first presented in Ref. [13] (see also 
Ref. [18]). 

We begin by using the fact that, on a foliage of space- 
time corresponding to spatially flat hypersurfaces, the 
curvature perturbation attributed to each of the universe 
components (labelled by the index i) is given by [25] 



-H 



(16) 



where pi and 6pi are, respectively, the energy density and 
its perturbation of the component in question. 

The total curvature perturbation C(i), which is also 
given by Eq. (16) with pi and Spi replaced, respectively, 
by the total energy density of the universe p = pi and 
its perturbation Sp, may be calculated as follows. Using 
the fact that 5p = J2i ^Pi ^^'^ continuity equation 
Pi = —3H(pi +Pi), where pi is the pressure of the i-th 
component of the universe, it is easy to find that 



Pi 



p+p 



Ci, 



(17) 



where p = '^iPi is the total pressure. Now, since in 
the curvaton scenario, all contributions to the curvature 
perturbation other than the curvaton's are negligible, we 
find that 



c = c. 



1 



Pa_ 
P 



(18) 



where ^ ~ 2 x 10~^ is the curvature perturbation ob- 
served by the cosmic microwave background explorer 
(COBE) [26], Wa and w are the curvaton and the overall 
barotropic parameters respectively (with Wi = Pi/pi) and 
Cct is the partial curvature perturbation of the curvaton. 
The right hand side (RHS) of this equation is evaluated 
at the time when the curvaton decays and this is indi- 
cated by the subscript "dec" . This decay occurs after the 
end of inflation in which case Eq. (18) gives 



(19) 



where fidec is the ratio of the curvaton energy density to 
the total energy density of the universe at the time of the 
decay of the curvaton: 



O = 

^ ^dec — 

p 



< 1. 



(20) 



From the bound [27] on the possible non-Gaussian com- 
ponent of the curvature perturbation from the recent 

cosmic microwave background radiation (CMBR) data 
obtained by the Wilkinson microwave anisotropy probe 
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(WMAP) satellite, one finds [25] that, at 95% confidence 
level (c.l.). 



10-2 < ridec < 1. 



(21) 



The partial curvature perturbation of the curvaton when 
the latter oscillates in a quadratic potential [cf. Eq. (15)] 
is given [24] by 



(a 



2 5a 

3 ~ 



a 



(22) 



where "osc" denotes the time when the curvaton oscilla- 
tions begin. 

In this paper, we assume that the Hubble parameter 
during inflation is comparable to the (tachyonic) masses 
that the radial fields which determine the value of the or- 
der parameter of the PNGB acquire after inflation. This 
means that the evolution of the curvaton's order parame- 
ter V ceases at (or soon after) the end of inflation. There- 
fore, at the end of inflation, v Vq and the mass of the 
curvaton assumes its vacuum value m^. Hence, in the 
following, we assume that the curvaton mass has already 
assumed its vacuum value before the onset of the curva- 
ton oscillations. Consequently, the curvaton oscillations 
begin when 



Ho 



1 



V3 



(23) 



Before the oscillations begin, the phase corresponding to 
the curvaton degree of freedom is overdamped and re- 
mains frozen. More precisely, this means that 



S9o 



(24) 



where the subscript star denotes the values of quantities 
at the time when the cosmological scales exit the horizon 
during inflation. 



a 

V 



(25) 



with 6 G (— TT, tt] and SO is its perturbation. Hence, for 
the curvaton partial perturbation, we find 



69 


_ se 


5a 









5_a_ 
a 

Now, for the perturbation of the curvaton, we have 

2^' 



(5(7* 



(26) 



(27) 



We assume that the order parameter of the PNGB during 
inflation is smaller compared to its value in the vacuum 
by a factor 



e = — < 1. 
Combining Eqs. (24)-(28), we find that 

2'K£ 



(28) 



(29) 




FIG. 1: Schematic representation of the amplification of the 
PNGB curvaton perturbation when the order parameter v in- 
creases from the value it has when the cosmological scales 
exit the horizon «, = evo to its vacuum value vq. V is the po- 
tential of a fiducial complex field <j} parametrizing the growth 
of the order parameter such that v/\^ — The phase 

6 oi (j) corresponds to the PNGB degree of freedom a (i.e. 
6 = (t/v). The perturbation in cr at horizon crossing has 
amplitude (5(t, ~ H^, which corresponds to a phase pertur- 
bation of magnitude 56 = 6a»/Vf As the order parameter 
grows, 56 remains constant (the phase perturbation is frozen 
on superhorizon scales) but the amplitude of the curvaton 
perturbation is increased up to 5<j ~ £~^J/*. 



which means that after the end of inflation, when v as- 
sumes its vacuum value, the curvaton perturbation is am- 
plified by a factor e^^ (see Fig. 1). From Eqs. (19) and 
(22), we have 



c 



■ 5ao 



Using Eq. (29), we can recast the above as 

We may obtain a lower bound on s as follows: 

6a^, iJ, 
< 1 ^ e > £m 



(30) 



(31) 



(32) 



< 



cr* 2ttvo 
where we have used Eqs. (26) and (29) and that ctq 

B. The bounds on the inflationary scale 



As is shown in Ref. [13], in the case when the curvaton 
oscillations begin after the order parameter of the PNGB 
has attained its vacuum value, we have 



dec 



min{mc,,rinf} 
max{iJdom, Po-} 



-ffBBN 



(M^Ti^^)-^, (33) 
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or equivalently (using V*^'^ ~ -\/3^/^iI*Mp) 



min{mo-,rinf} 



^(^^^^§^^Y {Ml,T^..)K (34) 

V -"BBN / 

where Finf and T„ are the decay rates of the inflaton and 
the curvaton fields respectively, -ffdom is the Hubble pa- 
rameter at the time when the curvaton energy density 
dominates the universe (if the curvaton does not decay 
earlier), and i/eBN = (tt / 2,) ^ g^^^ / IQ (TIbn/^^p) and 
^BBN ~ 1 MeV are, respectively, the Hubble parameter 
and the cosmic temperature at the time of big bang nucle- 
osynthesis (BBN) with .gBBN = 10.75 being the effective 
number of rclativistic degrees of freedom at that time. 

Now, we require that the curvaton field decays before 
BBN, i.e. To- > JIbbn- We also have Finf < -ff*. Hence, 
Eqs. (33) and (34) provide the following bounds 

> l^de! W^(^p2^bbn)*~ (^) ' X 10^ GeV 

(35) 

and 



bounds in Eqs. (35) and (36) with those in Eqs. (39) and 
(40) respectively, we find that the first set of bounds is 
more stringent if 



(41) 



Thus, for £ <C 1, the second set of bounds is typically less 
stringent than the first one. 

From Eqs. (32) and (36) and after a little algebra, it is 
easy to get 



> (^\^ 10-13 
^ if, >f^V 10-44 G^v^ 



(42) 



which means that, in principle, the larger vq is the 
smaller v}^^ can become. 



C. Scale invariance requirement 



1 / e 



dec 



X 10^2 Qgv. 



Furthermore, we also note that, generically, 
F >^ 



(36) 



(37) 



where the equality corresponds to gravitational decay. 
The above can be shown [13] to imply that 

if* > fidecM'Mp (I;) (max|l,^|y, (38) 



which results in the bounds 



H. > SiiU^EO'Mp (9^ 



and 



^ ) X 10^4 Qev. (40) 



These bounds may be relaxed if e is small enough. For 
a PNGB curvaton in particular, we may have mo- <IC i?* , 

which also relaxes the bounds in Eqs. (39) and (40). How- 
ever, in our case, rua ~ if* (see below). Comparing the 



The evolution of the order parameter v{t) during infla- 
tion is subject to an important constraint which has to 
do with preserving the scale invariance of the spectrum 
of the curvature perturbations. 

The amplitude of the curvature perturbation is deter- 
mined by the magnitude of the perturbation in the cur- 
vaton field, which, in this scenario, apart from the scale 
of Hif is also determined by the amplification factor . 
The latter is determined by the value of the order pa- 
rameter t;* when the curvaton quantum fluctuations exit 
the horizon during inflation. A strong variation of v{t) 
at that time results in a strong dependence of e{k) on 
the comoving momentum scale fc, which would reflect it- 
self on the perturbation spectrum threatening significant 
departure from scale invariance. 

In Ref. [13], it was shown that, in order for this to be 
avoided, the rate of change of the order parameter must 
be constrained as 



(43) 



From the above, it is evident that, in order not to violate 
the observational constraints regarding the scale invari- 
ance of the curvature perturbation spectrum, the order 
parameter must either remain constant or, at most, have 
a very slow variation when the cosmological scales exit the 
horizon [28] . However, this cannot remain so indefinitely 
because we need f o ^ to have substantial amplifica- 
tion of the perturbation (i.e. e <C 1). Consequently, v 
has to increase dramatically at some point after the exit 
of the cosmological scales from the inflationary horizon. 
This requirement is crucial for model building. 
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In this paper, wc will show that the evolution of v 
begins at a phase transition during inflation, as presented 
in Ref. [18]. Initially, the growth of v is very slow, but 
later, near the end of inflation, v grows substantially until 
it reaches its vacuum value vq. 



IV. CAN WE CONSTRUCT A PQ MODEL 
WITH A PNGB CURVATON? 



We will now address the question whether we can con- 
struct a PQ model [19] which, in addition to the standard 
axion, contains another axion-like field with the right 
properties to play the role of a PNGB curvaton as dis- 
cussed in the previous section. It is well known that, 
in the SUGRA extension of the minimal supersymmet- 
ric standard model (MSSM), there exist certain D- and 
F-flat directions in field space which can generate inter- 
mediate scales 



symmetry: 



(44) 



where n is a positive integer. It seems natural to try to 
identify Mi with the symmetry breaking scale /a of the 
PQ symmetry U(1)pq, such that a /i term is generated 
with fi ~ f2^^/M^ - m3/2 [29]. This would simultane- 
ously resolve the strong CP and ji problems of MSSM. 

The resolution of the /i problem forces us to consider 
non-renormalizable superpotential terms such as 



(45) 



where A is a dimensionless paramciter, P is a SM singlet 
superfield and h\, /12 are the electroweak Higgs doublets. 
The fact that the PQ symmetry carries QCD anomalies 
implies that the combination must have a non-zero 
PQ charge as one can easily deduce from the Yukawa 
couplings of the quarks. The field P must then necessar- 
ily carry a non-zero PQ charge. So, if it acquires a VEV 
of order Mi, the PQ symmetry breaks spontaneously and 
a n term of the right magnitude is generated via the su- 
perpotential term in Eq. (45). 

However, the field P has no self-couplings due to its 
non-zero PQ charge. Moreover, its couplings to the 
MSSM fields involve at least two of them since there are 
no SM singlets in MSSM. As a consequence, before the 
electroweak symmetry breaking, P has a fiat potential. 
To lift the flatness of its potential and generate an inter- 
mediate VEV for P of the order of Mi in Eq. (44), we 
must introduce [30 32] a second SM singlet superfield Q 
with non-zero PQ charge having a coupling of the type 



^P 



.n+3 — fc 



Q^/M, 



(46) 



where ^ is a dimensionless coupling constant and k is a 
positive integer smaller than n -|- 3. The superpotential 
term in Eq. (46) determines the PQ charge of Q. 

After soft SUSY breaking, wc obtain the following 
scalar potential for the spontaneous breaking of the PQ 



V = [(n + 3-fc)2|Q|2 + fc2|p|2] 



pn+2—kQk—l 



M^ 



+mj,\P\' + ml\Q\' 



+ 



A C ^ + h.c. 



(47) 



ip, m,Q ~ "^3/2 and can have either sign, while A 
is a complex parameter with magnitude of order 7713/2- 
For large enough this potential possesses non- 

trivial (local) minima at 



\P\,\Q\^{m3/2M^y' 



(48) 



even if mp, tuq are positive since the last term in the 
RHS of Eq. (47) can be adequately negative. We see 
that here the PQ scale /a is generated dynamically and 
is not inserted by hand as in the PQ schemes with renor- 
malizable interactions. 

In order to implement our scenario, we need a valley of 
local minima of the potential (with respect to the direc- 
tion perpendicular to the valley) which has negative incli- 
nation. Along this valley, the fields |P| and \Q\ must take 
values which arc much smaller than their vacuum values. 
If the system happens to slowly roll down this valley dur- 
ing the relevant part of inflation, the order parameter v of 
the PNGB remains, during inflation, much smaller than 
its vacuum value vq and our amplification mechanism 
for the curvaton perturbations may work. This can be 
achieved only if one of the masses-squared mp, TOq is 
negative. Let us assume, for dcfiniteness, that mp < 
and niQ > 0. Note, however, that the following discus- 
sion appHes equally well to the opposite case too, where 
TOp > and TOq < 0. In the case under consideration, 
the scalar potential is [32] unbounded below on the P 
axis (i.e. for (5 = 0) unless fc = 1 since, for A; > 1, all the 
terms in the RHS of Eq. (47) vanish on the P axis except 
the negative mass term of P. So, we restrict ourselves to 
the case k = 1. 

During inflation, m%, ttiq acquire [2, 14, 33] SUGRA 
corrections of order which arc assumed to be positive. 
Also, A receives SUGRA corrections of order H. As al- 
ready mentioned, the Hubble parameter during inflation 
is, in our case, of order So, in the initial stages of 

inflation, the effective mass-squared of P after SUGRA 
corrections, which we will call mp, can be positive. In 
this case, the origin in field space becomes a local mini- 
mum and the system may be initially trapped there. As 
H gradually decreases during infiation. Trip becomes neg- 
ative and the system starts slowly rolling down in the P 
direction. For non-zero P, the last term in the RHS of 
Eq. (47) yields a linear term in Q (recall that k = 1) and, 
thus, Q is also shifted from zero. More precisely, we get 



\Q\ 



|p|n+l 



W3/2 



M^ 



\P\ « \P\. 



(49) 
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The last term in the RHS of Eq. (47) then yields 



V2\P\ V2\Q\ 



(50) 



where (j)p and (^q arc canonically normalized real fields 
corresponding to the phases of P and Q respectively and 
A is the effective A after SUGRA corrections. Note that, 
in deriving Eq. (50), the fields P and Q were appropri- 
ately rephased so that the product A^ is negative. This 
is a convenient choice since, in this case, the expression 
in the above equation is minimized when the argument 
of the cosine vanishes. The orthogonal axion direction, 
which we would like to use as a PNGB curvaton, corre- 
sponds to the canonically normalized real field 



{n + 2)\Q\^P + \P\<j>Q 
V(n + 2)2|Q|2 + |P|2- 



(51) 



Its mass-squared can be evaluated from the term in 
Eq. (50), which, in view of Eq. (49), yields 



M^\Q\ ""3/2- 



(52) 



So, the mass of the orthogonal axion during inflation is of 
order 171^/2, which is comparable to the inflationary Hub- 
ble parameter. Consequently, this field docs not qualify 
as a curvaton, because it cannot obtain a superhorizon 
spectrum of perturbations. In conclusion, we have seen 
that our scenario cannot be realized within extensions of 
the MSSM with a PQ symmetry which contain only two 
SM singlet superfields. 

The addition of a third SM singlet superfield S, how- 
ever, can drastically change the situation allowing the 
implementation of our mechanism. We could keep the 
masses-squared of P and Q positive and include a su- 
perpotential term of the type in Eq. (46) with any value 
of the integer k between unity and n + 2. In this case, 
as already mentioned, the potential for P and Q (with 
S = 0) can possess non-trivial minima at \P\ and \Q\ 
given by Eq. (48). Now, we introduce an extra superpo- 
tential term of the type 



^^p^+3-p-qQPS1/M^, 



(53) 



where p, q are non-negative integers with p + q < n + 3 
and q> 3. This term determines the PQ charge of S. Of 
course, we should keep in mind that all possible terms 
involving P, Q and S that satisfy the global symmetries 
(including R symmetries) of the terms in Eqs. (46) and 
(53) should be present in the supcrpotential. We assume 
that the term in Eq. (53) is the term of this type with 
the smallest power of S. We take the mass-squared of S 
negative. Then, for small values of S, we obtain a valley 
of minima with negative inclination at almost constant 
values of \P\ and \Q\ given by Eq. (48). 

If the system slowly rolls down this valley during in- 
flation, the A-tevm corresponding to the coupling in 



Eq. (46) generates a mass term for the canonically nor- 
malized fleld 



{n + 3-k)\Q\(Pp + k\P\(j)Q 
V(n-F3-A;)2|Q|2 + fc2|p|2 



(54) 



with mass-squared of order TOg^j- This is easily shown by 
repeating the argument which led to Eqs. (51) and (52). 
The j4-term corresponding to the coupling in Eq. (53) 
yields 



- 2|i||^,| 



|p|n-K3-p-g|Q|p|5^|, 



X cos 



{n + 3-p- q)(j)p p(j)Q 



x/2|P| 



-I- 



\/2|Q| ' V2\S\ 



+ 



,(55) 



where 05 is a canonically normalized real scalar field cor- 
responding to the phase of S and was taken negative 
by repliasing S. This generates a "mass term" for the 
canonically normalized field 

{n^3-p-q)\Q\\S\cl>p+p\P\\S\^Q + q\P\\Q\d,s 

V'(n + 3-p-g)2|g|2|5|2+p2|p|2|5|2+^2|p|2|Q|2 

(56) 

with "mass-squared" 



q^\A\M\Pr+^-P-'^\Q\P\Sr 
MS 



~ m 



\S\ 



3/2 I IP 



q-2 



(57) 



for 1 5*1 ^ \P\ ^ \Q\ (we use quotation marks to indicate 
that the "mass term" and "mass-squared" do not neces- 
sarily correspond to a mass eigenstate). We see that, for 
g > 3, this "mass-squared" is suppressed relative to rn^^^- 
Supcrpotential terms with higher powers of S obviously 
give even more suppressed "masses-squared" . So, the lin- 
ear combination of 0p, (f>Q and (j)s which is orthogonal 
to the combination in Eq. (54) and the axion direction 
has mass-squared much smaller than m'^^^ during infia- 
tion and can be used as PNGB curvaton. It is though 
important to make sure that at least one of the three 
terms S^+^/Mg, S'^+^P/M^ and S'^+^Q/M^ is allowed 
in the superpotential, since otherwise the potential will 
be unbounded below on the S axis. In the next section, 
we will present a concrete class of models of this category. 



PQ MODELS WITH AN AXION-LIKE 
CURVATON 



We consider a class of simple extensions of MSSM 
which are based on the SM gauge group Gsm, but also 
possess two continuous global U(l) symmetries, namely 
a PQ symmetry U(1)pq and a R symmetry U(1)r., and 
a discrete Z2' symmetry. In addition to the usual MSSM 
left-handed superfields hi, /12 (Higgs SU(2)l doublets), 
li (SU(2)l doublet leptons), (SU(2)l singlet charged 
leptons), qi (SU(2)l doublet quarks), and u^, rf? (SU(2)l 
singlet anti-quarks) with i = 1,2,3 being the family in- 
dex, the models also contain the SM singlet left-handed 
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supcrficlds P, Q and S. The charges of the superfields 
under U(1)pq and U(1)r arc 

PQ: P(-2), g(2), 5(0), /ii, /i2(n + 1), 
R: P(^), Qi^), S{^^), h,M{Q) (58) 

with the "matter" (quark and lepton) superfields having 
PQ = -(n + l)/2, R = (n + l)(n + 3)/4. The integer n 
is taken to be of the form 

n = 4Z + 1, (59) 

where ? = 0, 1, 2, ... is any non-ncgativc integer providing 
a numbering of the models in this class (see Sec. VI B), 
and the charges are normalized so that they take their 
absolutely smallest possible integer values. Finally, under 
the Z2 symmetry, P changes sign. 

The most general superpotential compatible with these 
symmetries is 

W = yeij{lihi)e^j + yuij{qih2)u'=j + ydij{qihi)d^j 
+XP''+\hih2)/M^ 

(n+3)/4 

+ ^ Afe5"+3-^'=(P0)2'=/M^, (60) 

fc=0 

where y^j, Vuij, Vdij are the usual Yukawa coupling 
constants. A, Xk arc complex dimensionless parameters, 
{XY) indicates the SU(2)l invariant product CabXaYi, 
with e denoting the 2x2 antisymmetric matrix with 
£12 = 1, and summation over the family indices is im- 
plied. The R charge of W is (n-|-l)(n-|-3)/2. Baryon 
number is automatically conserved to all orders in per- 
turbation theory as a consequence of the R symmetry. 
The reason is [34] that the R charge of any combination 
of three color triplet or anti-triplet superfields exceeds the 
R charge of W and there are no superfields with negative 
R charge to compensate. Note that the Z2 subgroup of 
U(l)pQ coincides with the discrete matter parity sym- 
metry (denoted by Z^^), which changes the sign of all 
matter superfields. It is obvious that the superpoten- 
tial in Eq. (60) also conserves lepton number, which is a 
consequence of both the R and the PQ symmetry. 

Note that the superpotential in Eq. (60) is of the same 
type as the superpotential which has been considered in 
Ref. [10]. The main difference is that here the discrete 
Z„_|_3 symmetry of Ref. [10] is replaced by a much more 
powerful continuous U(l) R symmetry and, also, that an 
extra discrete Z2 symmetry is added. Moreover, in con- 
trast to Ref. [10], we include here all the superpotential 
terms which are compatible with the symmetries of the 
model. In this sense, the superpotential in Eq. (60) is 
a completely natural superpotential. It should be em- 
phasized that continuous global symmetries such as the 
U(1)r or the U(1)pq symmetry used here, rather than 
being imposed, can arise in a natural manner from an 
underlying superstring theory. Indeed, as it was pointed 



out in Ref. [35], discrete symmetries (including R sym- 
metries) that typically arise after compactification could 
effectively behave as if they are continuous. 

To sec that the above superpotential has the most gen- 
eral form allowed by the symmetries, observe that, due to 
matter parity, any term in W must contain an even num- 
ber of matter superfields. Actually, we may have either 
no or two matter fields since higher combinations carry R 
charge larger than that of W. The possible combinations 
of two matter fields are of the type II, e^e'^, Ze°, qu" and 
qd'^ (from color conservation) and have the R charge of 
W. So, we can multiply them only by superfields of zero 
R charge, i.e. hi, h^, and Q in the case n= 1. However, 
multiplying II, e'^e'^ just by Q's, we cannot compensate 
their non-zero weak hyperchargc. We need to multiply 
them at least once by hi or Indeed, the weak hy- 
perchargc of II (e'^e'^) is compensated if we take the com- 
bination llh^ha (e'^e'^hihihihi) whose PQ charge is n-|-l 
(3(n+l)), which could only be cancelled by including P's, 
the only superfields with negative PQ charge. This is 
though not allowed by R symmetry. So, the only mat- 
ter field bilinears which are allowed are qu", qd" and le^, 
which conserve lepton number. Actually, to cancel all 
the SM quantum numbers and the PQ charge, we must 
take the combinations {qh2)u'', {qhi)d'^ and {lhi)e'^. No 
further superfields can be included in these combinations 
since, by U(1)pq, we could only include the combination 
PQ, which has though positive R charge. In conclusion, 
we see that the only superpotential terms involving mat- 
ter superfields which can be present are the usual Yukawa 
couplings. 

We still have to consider superpotential terms with no 
matter superfields. If such terms involve hi, h2, these 
superfield must enter through the combination {hih2) 
which is neutral under both SU(2)l and U(l)y. To cancel 
the PQ charge, we must then take P"'~^^ {hih2) ■ To retain 
the U(1)r symmetry, we could multiply this only by the 
combination PQ, which has PQ = 0. However, the R 
charge of P"+^(/iift.2) is already equal to the R charge of 
W and PQ has positive R charge. So, P"+^ (/iift,2) is the 
only allowed term which involves Higgs but no matter 
superfields. We are left to consider combinations which 
involve only SM singlet superfields. The PQ symmetry 
allows only the combinations Sp{PQY, where p and q 
are non-negative integers. The R charge of PQ (S) is 
n + 1 ((n + l)/2), which implies that the term 5"+^ can 
be present in the superpotential. Note that the combina- 
tion PQ has the same R charge as S^. Consequently, the 
terms S'"+i(PQ), ^"-^(PQ)^,..., (Pg)("+3)/2 are also 
allowed by U(1)r. The symmetry, however, forbids 
all the odd powers of PQ and we arrive at the terms in 
the sum in the RHS of Eq. (60). 

The PQ symmetry is anomalous (as it should). In 
particular, one can show that the QCD instantons break 
it explicitly to its Zj^ subgroup with J\f = 6{n + 1). 
This subgroup contains the Z^^ (matter parity), which 
thus always remains unbroken by instanton effects. On 
the contrary, the R symmetry is non-anomalous since 
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the fermionic components of all the color triplet or anti- 
triplet superfields (actually, all the matter fermions) have 
zero R charge. The soft SUSY-breaking terms (especially, 
the ^-type terms), however, break U(1)r to its Zm. sub- 
group, where = (n + l)(n + 3)/2. 



VI. THE SCALAR POTENTIAL 

The part of the superpotential in Eq. (60) which is 
relevant for the PQ breaking is the sum in the RHS of 
this equation. This sum contains at least two terms. 
For n = 1, in particular, it consists of just two terms. 
The resulting scalar potential after including soft SUSY- 
breaking effects is 



where 



V = \Fi 



(n+3)/4 



(61) 



S'"+3-4'=(pg)2fe-iQ 



fe=l 



FQ, (62) 



(n+3)/4 cn+3-4fe/priA2fc-l p 

Fq= 2fcA.^^ ^I^^^^FP (63) 



and 



fc=i 



(n-l)/4 



Fs= {n + i-Ak)Xk- 



™+2-4fe 



{PQ) 



2k 



are the F-terms, and 

14oft = ml\P\'' +ml\Q\'' +m%\S\^ 

(n+3)/4 cn+3-4fe/'pr)^2fc 
fc=0 



(64) 



(65) 



the soft SUSY-breaking terms. Here, the soft SUSY- 
breaking masses-squared rup, rag and m| are of the 
order of the ^^^2 ^^'^ can have cither sign. Also, for 
simplicity, we assumed universal soft SUSY-breaking A- 
terms with the magnitude of the complex parameter A 
being of the order of m3/2. Note that the sums in the 
RHS of Eqs. (62)-(64) contain at least one term. Actu- 
ally, these sums consists of just one term for n = 1. On 
the other hand, the sum in Eq. (65) contains at least two 
terms with the minimum number of terms corresponding 
to n = 1. 

As |P|, l^l 00, the potential is generally dom- 
inated by the positive F-terms and thus V is bounded 
below independently of the sign of the soft mass terms. 
However, on the P or Q axis, i.e. for Q = 5 = or 
P = 5 = respectively, the F-terms vanish identically 
together with the A-terms and the potential is given by 
just the mass term of P or Q respectively. So, to have V 



bounded below on these two axes too, we must restrict 
the masses-squared of P and Q to be positive. For sim- 
plicity, we will take these two soft masses-squared to be 
equal, i.e. we will put 



Top = to: 



Q 



TO . 



(66) 



The potential on the S axis (i.e. for P = Q = 0) is 



y = (n+3)2|Ao|' 



2(n+2) 



Ml" 



-+ml\S\ 



S 



n+3 



Mi 



h.c. 



(67) 

with the first term in the RHS of this equation originating 

from the F-term Pg. We see that, for \S\ — > oo, the 
potential is dominated by the positive F-term and thus 
V on the S axis is bounded below no matter what the 
sign of the mass-squared of S is. So, this sign can be 
chosen at will. For reasons which will become clear later, 
however, we take 



to| < 0. 



(68) 



Therefore, the origin in field space (P = Q = 5 = 0) is 

a saddle point of the potential with positive curvature in 
the P and Q directions and negative in the S direction. 

Using Eqs. (61)-(65), one can readily show that, for 
TO.p, nig > 0, the potential V in Eq. (61) has a valley of 
local minima with respect to |P| and \Q\ which lies on the 
S axis (i.e. at P = Q = 0). The potential right on the 
bottom line of this valley can be found from Eq. (67) by 
choosing the phase 9s of S so that the sum of the terms 
in the parentheses in the RHS of Eq. (67) is minimized. 
In this case, this equation takes the form 



|5'|2(n+2) 



V = {n + sriXol'^^^^ + mi\S\' - 2\A\\\o 



^s\ 



n+3 



which has a minimum at 



1^1 



n+l 



Mi 



\A\ + ^|A|2 -4(n-h2)TO^ 
2(n + 2)(n + 3)|Ao| 



(69) 



(70) 



where IS*] ^ (m3/2Mp)~. So, as \S\ increases from 
zero, the depth of the valley increases (i.e. the valley has 
initially negative inclination) until l^j reaches the value 
in Eq. (70), where the maximal depth is achieved. As \S\ 
increases further, the bottom line of the valley rises and, 
for \S\ oo, tends to infinity. We will call this valley, 
which starts from the trivial saddle point at the origin, 
the trivial valley. 



A. The non-trivial minima of the potential 

The extrema of the full potential V in Eq. (61) are 
given by the equations 

dV _ dFp OFq , p*5Ps 2p* , ,p _ f) 



(71) 
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dQ ^ dQ ^ dQ 



and 

dV _ dFp OFq 
95 -^^^+^^3^ 



F*c:^+m^Q* +AFq = Q 



dQ 



(72) 



F*s^+m%S*+AFs = 0. 

(73) 

Multiplying Eqs. (71) and (72) by P and Q respectively, 
subtracting and using Eqs. (62)- (64), we obtain 



(|F|2 + m2)(|P|2_|0|2)=0, 



(74) 



which implies that |P| = \Q\. So, the complex fields P 
and Q have exactly the same magnitude in any extremum 
of the full potential V. This is actually true also in the 
cxtrcma of V with respect to P and Q only for fixed S 
since wc have not used Eq. (73). In the trivial minimum 
in Eq. (70) which lies on the trivial valley, \P\ = \Q\ = 0. 
However, the potential V possesses non-trivial minima 
too, where 



|Q|,|5|~(m3/2A/f^ 



(75) 



It is important to note that the continuous global sym- 
metry of the model is not exactly U(1)pq x U(1)r. The 
reason is that there are elements of U(1)r which are in- 
distinguishable from elements of U(1)pq since they have 
the same action on all the superfields of the model. Actu- 
ally, the only elements of U(1)r which can, in principle, 
be identified with elements of U(1)pq are the ones be- 
longing to its Zm subgroup with M. = (n -|- l)(n -|- 3)/2 
generated by the element 



e'^ e U(1)r. 



(76) 



where Zn+s is generated by 



e'TTfs e U(l)r 



(80) 



So, the global symmetry of V takes the simple form 
U(l)pQ X Z„+3 X Z^. 

In any non-trivial minimum (actually, for any fixed 
non-zero values of P, Q and 5), this symmetry is spon- 
taneously broken to the Z™^ subgroup of U(1)pq. In 
case there was, after inflation, a phase transition from 
the origin in field space (P = Q = 5 = 0) to a non-trivial 
minimum, we would encounter copious production of ax- 
ionic strings [36] as well as domain walls. As we perform 
a full rotation around such a string, the phase of P or 
Q changes, respectively, by —2n or 27r. There are two 
types of walls separating vacua which are related either 
by the group element in Eq. (80) or the generator of Z2' ■ 
Note that n + 3 walls of the former type can terminate 
together on the same line. 

As explained in the previous section, after the onset 
of instantons at the QCD transition, U(1)pq is explicitly 
broken to Zj^ with M = 6{n + 1) which is generated by 
the element 



e'ei^r+T) g U(l)pQ. 



(81) 



So, after instantons, the axionic strings become [36] 
boundaries of 3 (n -I- 1) axionic walls [37] which separate 
vacua related by the group clement in Eq. (81). We see 
that, if, after inflation, a transition from the origin in 
field space to a non-trivial minimum takes place, a rich 
system of domain walls is produced leading to a cosmo- 
logical catastrophe. So, it is clear that such a transition 
should be avoided in our model. 



This is so because this Z^ is the maximal ordinary (i.e. 
non-R) symmetry group contained in U(1)h. One can 
show that the element 



B. The shifted valley of minima 



e u(i)r 



(77) 



has the same action on all the superfields as the element 
e'^ e U(l)pQ (78) 

and, thus, the ^(„+i)/2 subgroup of Z^ which is gen- 
erated by it is identical with the Z(^n+i)/2 subgroup of 
U(l)pQ generated by the element in Eq. (78). So, the 
continuous global symmetry of the model is actually 

U(1)PQ X (U(1)r/Z(„+i)/2). 

As already explained, U(1)r is a symmetry of the 
model only in the limit of exact SUSY. In the scalar po- 
tential V which includes the soft SUSY-breaking terms 
too, it is explicitly broken to its discrete subgroup Zj^. 
Thus, the group of global symmetries of V (except ) 
is U(l)pQ X {ZM/Z^n+i)/2)- The Zm symmetry can be 
factorized as follows: 



(79) 



For l^l < jPj - jgj, the F-terms in Eqs. (62)-(64) can 
be approximated as follows: 



n + 3^ {PQ) — Q 

n-1 S\PQffQ 
+ ^^^ + •••' (^2^ 



n + 3^ (PQ)^P 

An+3 ■ 



Fq = —^Xn±3 

n — 1 



S^{PQ)^P 



+ ..., 



(83) 



(84) 
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while the expression in the brackets in the RHS of 

Eq. (65) takes tlie form 

-1 + 3 



A\r 



{PQY 



+ AX,_ 



S\PQ) — 



+ h.c. \+...., 



(85) 

where the eUipses in these equations represent terms of 
higher order in S. Note that the second term in the RHS 
of Eqs. (82) and (83) exists only for n > 5 (/ > 1), while 
the second term in Eq. (85) exists for all values of n in 
Eq. (59). Using these relations, we can expand, in this 
regime, the potential V in Eq. (61): 



V = 



(n + 3)2 



+m%\P\ 



,\PQY'+\\P? + \Q?) 



X cos 



+ m||g|2-2|A||A^ 
f 3 



{\p\\Q\y 



-{Op + Oq) 



+ ... 



/(o) -h . . . , (86) 

where 9p and 9q are the phases of the complex scalar 
fields P and Q respectively. Here we assumed, without 

loss of generality, that the product A\ri+3 is real and 

4 

negative, which can be readily achieved by redefining the 
phase of the product of fields PQ. 

The leading order part V(o) of the potential V in 
Eq. (86) (consisting of the explicitly displayed terms in 
the RHS of this equation) is minimized with respect to 
the phases Op and 9q for 

n + 3 

—^{6p+6q) = Q modulo 27r. (87) 

Under this restriction on the phases, the extrema of V(o) 
with respect to |P| and \Q\ are given by the conditions 



dV( 



(0) 



(n + 3)2 



d\P\ 



|A:i±3p[(n + 3)|P|2 + (n + l)|Qn 



mQ\ 



n+l 



Ml" 



■\Q\' 



2m^|P|-2|A|^^|A^| 



= (88) 



and 
d\Q\ 



(n + 3)2 



|A^r[(n+l)|P|2 + (n + 3)|g|2] 



|P|»+i|Q|" 



\P\ 



M2 



+ 2ml\Q\ 



= 0. 



2\A 



(n + 3) , 



(89) 



Multiplying Eqs. (88) and (89) by |P| and \Q\ respec- 
tively and subtracting, we obtain 

(n + 3)2 



-|A:^|^(|P|- 



+{^l\P? 



= 0, 



(90) 



which, for rrip = roq = > (see Eq. (66)), implies 
that |P| = \Q\. Substituting |P| for \Q\ in any of the 
Eqs. (88) and (89), we then obtain |P|=0 or 



(n + 3) 



p|2(n+l) 



^|A:4.|^(2n + 4)^^ + 
-2|^|l!^|A..3|^=0. 



2m^ 



(91) 



This equation has two real and positive solutions for 
\A\^ > 4(n-|-2)m^ 



which are given by 



|p|n+l 



x± 



|A[± v^|^P-4(n + 2)m2 
(n + 2)(n + 3)|An+3| 



(92) 



Obviously, x = and x = x+ correspond to (local) min- 
ima of V(o), while x = X- corresponds to a local maxi- 
mum (here x = |P|"+^/Mp). So, V(o) has a trivial min- 
imum at |P| = IQI = and a "shifted" minimum at 
|P| = \Q\ ~ (m3/2-Mp)~, which is the absolute mini- 
mum for \A\ > (n -h 2>)m. Note that the presence of the 
term An+s {PQ) ^ /Mp in the superpotential of Eq. (60) 

4 

is vital to the existence of the shifted minimum. In view 
of the Z2 symmetry, however, this superpotential term 
can only exist if (n -|- 3)/2 is an even positive integer, 
which implies the restriction in Eq. (59). 

The trivial minimum of V(o) cannot be consistent with 
our starting hypothesis that l^l <C |P| \Q\ and should, 
thus, be discarded. However, as we have already shown, 
it happens to exist as a minimum of the full potential V 
in Eq. (61) with respect to |P| and |Q| for all values of 
\S\ and constitutes the trivial valley of minima. 

To see how the shifted minimum of V(o) evolves as 
15*1 increases from zero, we consider the dominant S- 
dependent part of the potential V for l^j <C |P| ~ \Q\, 
which is given by 



AX, 



S\PQ)^ 
MS 



\n-3 



x(ipr+iQn ^™r +h.c. 



(93) 



where the first term in the brackets corresponds to the 
second term in Eq. (85), while the second term in the 
brackets originates from the interference of the two ex- 
plicitly displayed terms in the RHS of Eq. (82) plus the 
interference of the two explicitly displayed terms in the 
RHS of Eq. (83). So, the first term in the brackets ex- 
ists for all values of n, while the second only for n > 5. 
Note that V(i) contains the next-to-leading and the next- 
to-next-to-leading order parts of V in the expansion of 
Eq. (86), which are quadratic and quartic in S respec- 
tively. Actually, for |P| and at the shifted minimum 
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of V(o) J all the terms in V(o) are of the same order of mag- 
nitude, while the first (mass) term in the RHS of Eq. (93) 
is suppressed by (|S'|/|P|)^ and the terms in the brackets 
by (|5|/|P|)4.^ 

The minimization conditions for V with respect to 
\P\ and IQI coincide to leading order with Eqs. (88) 
and (89). The dominant corrections to these conditions 
for 15*1 <C |P| ~ IQI originate from the next- to-next-to- 
leading order terms in V^i) and are, thus, suppressed by 
(|5|/|P|)4. So, for l^l < |P| ~ \Q\, the shifted minimum 
of V(o) is also a minimum of V with respect to |P| and 
IQI at a practically S'-independent position. As a con- 
sequence, for small values of \S\, we obtain a "shifted" 
valley of minima of V at almost constant values of |P| 
and |Q|. This valley has obviously negative inclination 
for non-zero and small values of \S\, due to the nega- 
tive mass term of S. It starts from the "shifted" sad- 
dle point of V which lies at \S\ = and |P|, \Q\ equal 
to their values at the shifted minimum of V(o). Let us 
note, in passing, that a shifted valley of minima was first 
used in Ref. [38] as an inflationary trajectory in order 
to avoid the overproduction of doubly charged [39] mag- 
netic monopoles at the end of hybrid inflation in a SUSY 
Pati-Salam [40] GUT model. 



C. The PNGB curvaton 

The dominant S'-dependent part of V can be expressed 
in terms of the phases oi P, Q and S as follows: 



(1) 



ml\S\'^ -2\A\\Xr^ 



MS 



xcos(46»s + ^^^(6'p + 6'q) 



(n- l)(n-h3) 



X^\\X^\\S\\\P\' + \Q\') 



X cos {49s -2{ep + eQ)), (94) 

where we assumed that AA^^ is real and negative, which 

4 

can be readily arranged by redefining the phase of S. 
Note that all the other AXk's (except AXr^+s and AXn^ ) 

4 4 

remain in general complex since there is no extra field- 
rephasing freedom left. From the preceding discussion, 
we see that, on the shifted valley of minima of V, the 
second and third term in the RHS of Eq. (94) are of the 
same order of magnitude. For n = 1, the third term 
vanishes. Moreover, using Eq. (92), one can show that, 
for n > 5, the coeflicient of the cosine in the second term 
is always greater in absolute value than the coefficient of 
the cosine in the third term. Under these circumstances, 
V(i) is minimized with respect to the phases 6p, 6q, 6s 
of the fields by taking 



77 — 1 

^6s + ^{0p 



9q) = modulo 2-K. 



(95) 



This together with Eq. (87) implies that 

46is-2(6ip + 6»Q) = modulo 27r, 



(96) 



which maximizes the third term in the RHS of Eq. (94) . 

We can define the real canonically normalized fields 
corresponding to the phases of P, Q, S as follows: 



V2\p\ei 



= V2|Q|0Q, 0s = V2\s\es 



(97) 

On the shifted valley, the last term in the RHS of Eq. (86) 
generates a "mass-squared" for the real canonically nor- 
malized field (j)pQ = {(t)p + (j)Q)/V2 given by 



m 



PQ 



|^|(n + 3)^l^ ^^^^ (98) 



where = 
is of order 



|P|„- 

^3/2- 



"VA^p)+ is given by Eq. (92) and m.j,Q 
Also, the second and third term in 



the RHS of Eq. (94) generate on this valley 
squared" , respectively, for the combinations 

4V2|P|05 + (n-l)|^|0PQ 



masses- 



<f>S2 = 



v/32|P|2 + (n-l)2|5|2 
V2\P\<j>s - mpQ 



given, respectively, by 
m|i = 16|A||An^|a;+ 



(99) 



\S\' 



1 + 



32 



™S2 = -8(n- 



1 + 



l)(n- 

\S\' 
2|QP 



3)|A:^||A:43|4 



\sl 



(100) 



which are suppressed by |S'p/|(5p relative to m^^^. Note 

that the field a = {(j)p — (t>Q)/V2 remains massless to all 
orders in perturbation theory since it corresponds to the 
axion. So, we obtain a 2-dimensional subspace of mas- 
sive fields spanned by (f)pQ and (ps which are orthogonal 
to the massless axion direction. The quadratic form ob- 
tained by summing the three "mass terms" correspond- 
ing to the "masses-squared" in Eqs. (98) and (100) must 
then be diagonalized to find the two mass eigenstates and 
eigenvalues. This task is particularly simple in the limit 
|5| <C |P|. To lowest order in \S\/\P\, the eigenstates 
coincide with and (j)s, which we will also call a as 
it will be our PNGB curvaton (see below). The masses- 
squared of these fields on the shifted valley are equal to 
rripQ in Eq. (98) and 



ml 



\s\ 



2 r 



P 



in + 5)\A\ 



-(n- l)v'|A|2 -4(n-F2)m2 (101) 



respectively. Obviously, ifi^ is positive and suppressed 
by |S'|V|Pp relative to mj^Q, 



|P|2 relative to nipn, which is of order rn^^^- 
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the field 05 = (7 is a light PNGB when the system rolls 
down the shifted valley of the potential with \S\ <C \P\, 
while (j)pQ is a massive field. Inclusion of higher order 
corrections in the potential V along the shifted valley 
does not change this situation. 

The discrete Z2' symmetry of the model, under which 
P changes sign, is very important for the PNGB nature 
of (ps- Without this symmetry, the next-to- leading order 
term in the RHS of Eqs. (82) and (83) would be propor- 
tional to S"^ rather than S'^, the leading order term in 
the RHS of Eq. (84) would be linear in 5 rather than 
cubic, and the next-to-leading order term in Eq. (84) 
would contain 5^ instead of S"^. As a consequence, 
the terms in the dominant /S-dependent part of V for 
15*1 <C |P| ^ IQI which depend on the phases of the fields 
would be quadratic in S (compare with Eqs. (93) and 
(94)) and, thus, the mass of on the shifted valley 
would be of order 7723/2- 



D. Cosmological evolution 

SUGRA corrections [2, 14, 33] during inflation will add 
to ^ a term proportional to the Hubble parameter during 
inflation, which is of order 7713/2 in our case. To simplify 
the discussion, we take these corrections to be universal. 
Also, the masses-squared m%, tuq and rrtg will acquire 
corrections proportional to H^. We assume that these 
corrections are positive and, for simplicity, we also take 
them universal at least for mp and ttiq. So, nothing 
changes in the above discussion and formulae after in- 
cluding the SUGRA corrections during inflation except 
that A, and rrig must now be replaced by their effec- 
tive values 



A + caH, 



rn? = rr? -|- cpqH^ ^ m| = m| + csH^ 

(102) 

respectively. Here, ca is a complex parameter of order 
unity, while cpq and cs are real and positive parameters 
again of order unity. 

The effective parameters A and ?7i are of the same or- 
der of magnitude as A and m, i.e. they are of order 
7773/2. The mass-squared of S, which is taken negative, 
receives positive corrections from SUGRA, which are of 
the same order of magnitude as 7n|. We can arrange 
the parameters so that the effective mass-squared of S 
is positive during the initial stages of inflation. In this 
case, the shifted saddle point of V becomes a local min- 
imum of the effective potential and the system may be 
initially trapped in this minimum during inflation. As 
H decreases gradually during inflation, the SUGRA cor- 
rections become smaller and, at some moment of time, 
this minimum may turn into a saddle point. The system 
then slowly rolls down the shifted valley which has a very 
small slope given by the small effective mass-squared of 
S [41]. During this slow roll, (t)s is an effectively massless 
PNGB which can act as curvaton. 

After the end of inflation, the system keeps rolling 



down the shifted valley and eventually ends up in damped 
oscillations about a non-trivial minimum of V where all 
the fields P, Q and S acquire non-zero values of order 
(7ft3/2Mp)"+T and the PQ symmetry is broken. Note 
that, if, in the initial stages of inflation, the system hap- 
pened to be trapped in the trivial saddle point of the 
potential Va.tP = Q = S = Q (which, of course, is 
a local minimum of the effective potential in this case), 
it would later enter into the trivial valley along the S 
axis rather than the shifted one. So, it would end up 
in the minimum of Eq. (70), where P = Q = Q. This 
is obviously highly undesirable since, in this case, the 
PQ symmetry remains unbroken and no ji term is gen- 
erated. It is also important that, in our case, the fields 
P, Q and 5* have non-zero values during infiation after 
the time when the cosmological scales exit the horizon as 
they lie on the shifted valley. So, the global symmetry of 
the model is already broken to Z™'' during the relevant 
part of inflation. Consequently, neither domain walls nor 
axionic strings [36] are generated as the system settles 
in a non-trivial minimum of V . Also, no axionic walls 
[37] appear at the QCD transition since the spontaneous 
breaking of Zfj- to Z™"^ takes place before the relevant 
part of inflation. Therefore, no cosmological catastrophe 
is encountered. 



VII. CURVATON PHYSICS 
A. The curvaton potential 

From Eqs. (61), (84) and (94), we find that the domi- 
nant part of the scalar potential which is relevant to our 
curvaton candidate, in the case when \S\ <C |P| ~ \Q\, is 
given by 



= fhl\S\^-2\A\\\, 



l^l^(IJ^IIQI)^ 



X cos 



77 — 1 

^0s + ^^{6p + 6q) 



(77-l)(n + 3) 



\^\\\^\\s\Wp\^ + \q\^) 



x^^i||?^cos[4^s-2(^P + M] 
+16|A^| ^2n \S\ +••• , 



(103) 



where we included the SUGRA corrections during infia- 
tion and the ellipsis denotes terms of higher order in [S'l, 
which are, therefore, subdominant. In the above, accord- 
ing to Eq. (102), we have 



i^s = csH - \mg\ 



(104) 



where cs ^ +1 and |r7i|| ~ m? 



3/2- 



The potential Vcun is simplified by considering that the 
|P| and \Q\ fields have already assumed their minimum 
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value on the shifted vaUey as given by Eq. (92) with A 
and m replaced by A and m respectively. Furthermore, 
we can set the phases 9p and 0q equal to zero since we 
are only interested in the curvaton field direction, which 
practically corresponds to 9s for \S\ <C \P\. Then, in 
view also of Eq. (97), the above potential becomes 



+ 16\Xr_ 



K 2 COS 

K Ival 

|p|2(n-l) 
|2 l-^lval I CI 6 I 



4)3 



(105) 



where the value |P|vai of |P| on the shifted valley is given 
by l^'lvai V^p = a;+ (with x+ from Eq. (92), where A 
and m are replaced by A and m respectively) and 

(1 + Z) [(,; + .-,) - [a - l)Z\ 



K = 



Xr 



An+3 
4 

with Z being 



(n + 2)2(n + 3) 



Z= Wl-4(n + 2) 



m 



(106) 



(107) 



Note that |P|vai is practically constant when the cos- 
mological scales exit the inflationary horizon since the 
shifted valley is almost jS"! -independent for l^l <C 1^1 as 
shown in Sec. VI B and the Hubble parameter is very 
slowly varying (see Eq. (10)). 

From Eq. (107), it is evident that < Z < 1 and, there- 
fore, K is positive [42]. Hence, setting 



a :^ ipS and v 



1 



2V2 



(108) 



we obtain the curvaton potential for v <^ vq (the value of 
V in the vacuum) as 



V{a) ~ 64k|AP 



cos 



(109) 



Considering that the value of Vq is given by the value |S'|o 
of IS'I in the vacuum, for which |5|o ^ |-P|vah we find 



cos 



(9 



V{a) ~ 

where [cf. Eq. (44)] 

uo ~ Mi ~ (ms/sM^)^. 



(110) 



(111) 



Comparing the above with Eq. (15), we see that the mass 
of the curvaton is given by 



ma 



^'^1 (^) • 



(112) 



which agrees with Eq. (101). Therefore, when the cos- 
mological scales exit the horizon, we have 



rria ~ eTO3/2, 



where we used Eq. (28) and considered that \A\ ^ m3/2 
and K ~ 1. Hence, because, in the modular inflation 
model that we are considering, we have i?* ~ 1TI3/2, we 
find that, since e <C 1, u is, as required, effectively mass- 
less when the cosmological scales exit the horizon during 
inflation. 

The order parameter v for our curvaton fleld is de- 
termined by the value of \S\, for which the potential in 
Eq. (105), when taking ct — > 0, becomes 



V{\S\)c^fnl\S\ 



where 



-^Ival 



4(1 + Z) 



+ ■■■, (114) 



(115) 



and we have used Eq. (92). The minimum of the above 
potential occurs at 



Mi, 



(116) 



which, after the end of inflation, becomes equal to jiSjo 
as A and m are replaced by A and m respectively. 



B. The required e 

Let us now calculate the value of s required so that our 
curvaton scenario works. Firstly, we note that, in our 
case, the curvaton assumes a random value at the phase 
transition at which, during inflation, the system leaves 
the shifted saddle point of the potential and starts slowly 
rolling down the shifted valley (see Sec. VI D). This value 
typically is cr ^ u. After the end of inflation and before 
the onset of the oscillations, the phase 9 corresponding 
to the curvaton degree of freedom is overdamped and 
remains frozen. Hence, we expect that, at the onset of 
the oscillations, we have 



(117) 



where, typically, ~ ^5 ~ 1 and we took into account 
that the order parameter assumes its vacuum value very 
soon after the end of inflation. Combining Eqs. (31) and 
(117), we find 



77(9 V M] 



"^3/2 



(118) 



where we also used Eq. (Ill) and that iJ* ~ 7723/2- The 
e above is always larger than £min, where 



W3/2 
Mp 



(119) 



(113) which is derived from Eq. (32) with ~ ^3/2- 
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Let us now enforce the constraint in Eq. (35), which, 

for H^, ~ ra^ji, reads 



dec" 



(120) 



From Eqs. (118) and (120), it is easy to find that the 
above bound can be satisfied only if n is large enough: 



n > 



s + iog{nije) 
7-iog(fiL/^) 



(121) 



In view of Eq. (21), we see that, for ^ ~ 1, we have n > 1. 

An upper bound on n can be obtained by requiring 
that the curvaton decays before BBN. The interaction 
of a with ordinary particles is governed by the effective 
/i term in Eq. (45), which results [10] into the following 
decay rate of a into two Higgs particles: 



(122) 



Demanding that To- > Hbbn results in the bound 

~ 10-^ TeV > i/BBN ~ 10-2^ TeV 



=^ ma> 10— TeV, 



(123) 



where we used Eq. (111). The above requirement is al- 
ways satisfied if > 10 TeV. However, in the opposite 
case where rria < 10 TeV, we see that it yields an upper 
bound on n: 



n < 



9 + log(m^/TcV) 



(124) 



1 - log(m,/TcV) 
which, roughly, demands n < 9 for < 1 TeV. 

C. The reheating of the universe 

We will now proceed further by first discussing the 

reheating of the universe. This requires that we consider 
separately the cases when the curvaton decays before or 
after it dominates the universe. 

1. Curvaton decay before domination (i^dec "^^J 

During the radiation era and after the onset of curva- 
ton oscillations, for the curvaton energy density fraction, 
we have pa/p oc a{t) cx where a{t) is the scale fac- 

tor of the universe. Hence, in this case, wc find that 



ft. 



dec 



'^osc 



(125) 



where we have used Eq. (31) and 



Pa 

P 



Mv 



(126) 



which is derived from the fact that p, 
Pose - mlM^. Using Eq. (122) into Eq. (125) and also 
Eqs. (31) and (111), we obtain 



(127) 



Here, we have also used that Finf < i/* ~ 7713/2 ~ ma and 
Tinf ~ 5^m3/2 (128) 

with g being the dimensionless coupling constant of the 
inflaton to its decay products and the mass of the inflaton 
field s taken to be mg < i?* ~ ^13/2. 

In principle, g can be as low as mg/Mp if the inflaton 
decays gravitationally. However, since reheating has to 
occur before BBN, g has to lie in the range: 



10-" ^ 10!^ <g<l, 



(129) 



where we used the fact that the reheat temperature Treh, 
in this case, is 



Treh ~ VrinfMp ~ ^^ms/sMp. (130) 

Combining Eqs. (118) and (127), we find the relation 

71-2 

g 1 



02 V Mp 



(131) 



which results in 



n • 



m-\ogg + 2\og{^lJ6) 
lb + \ogg-2\og{nlJ6) 



(132) 



In view of Eqs. (21) and (129), we see that, for 6* ~ 1, the 
allowed range for n is 



2 < n < 44. 



(133) 



The lower bound in the above is tighter than the bound 
in Eq. (121). In fact, comparing Eqs. (121) and (132), it 
is easy to obtain the bound 



gm^<iQ'e, 



(134) 



which is easily satisfied provided that the angle 6 is not 
extremely small. 
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2. Curvaton decay after domination (Tidec w Ij 

In this case, the curvaton dominates the energy density 
of the universe when H = -ffdom, where -ffdom is given by 



dom 



-j^] m.m{m„,Tini}- (135) 



Now, using Eqs. (117), (122) and (128), it can be shown 
that the requirement < -ffdom results in the bound 



9 > 



1 / "l3/2\ 



which yields 



n > 



6»2 V Mp 

30 - log g - 2 log 6* 
15 + log 5 + 21og6'' 



(136) 



(137) 



This provides a lower bound on g for given n and 6, 
reminiscent of Eq. (132) with fidec ~ 1- For 6 1, the 
above bound implies 



n > 2. 



(138) 



This time the hot big bang begins after the decay of the 
curvaton, which suggests that the reheat temperature is 
now given by 



Treh ^ \/T„Mp 



Mp 



2 I. n+l J 



(139) 



It can be easily checked that the above is higher that 
Tbbn when n < 9, in agreement with Eq. (124). 



D. Avoiding axion overproduction 

There is a stringent upper bound on the PQ scale orig- 
inating from the requirement that the generated axions 
do not overdose the universe. The typical mass of the 
axion is [43] 



ma ^ 10 



10^2 GcV \ 

/a ) 



eV, 



(140) 



where, for the PQ scale (axion decay constant), we have 
/a ~ ^^0 [44]. The above mass implies that the onset of 
axion oscillations takes place when the energy density of 
the universe is 



10^ 



10^2 GcV 



GeV, (141) 



where the subscript "axosc" indicates the time at which 
axion oscillations begin. In contrast, the energy density 
of the universe when the curvaton decays is 



Pde 



10 



i^^) GeV, (142) 



Vo 



where wc have used Eq. (122) and also that nia 
Now, from Eq. (Ill), wc have 



which suggests that, in all cases, 

Pdec < 



ms/2. 
(143) 

(144) 



Hence, the axion oscillations always start before the cur- 
vaton decays. Thus, at the onset of the axion oscillations, 
the ratio of the axion energy density to the energy density 
of the universe is 



El 
P 



" \Mp) 



TO3/2 
Mp 



(145) 



where we have considered that the amplitude of the axion 
field at the onset of its oscillations is ~ ^a/a with Og, being 
the initial misalignment angle. 

Let us, firstly, investigate the case when the curvaton 
decays before it dominates the universe. In this case, 
the axion oscillates in a radiation background until the 
time teq of equal matter and radiation energy densities, 
denoted hereafter by the subscript "cq" . Since the energy 
density of the oscillating axion scales like pressureless 
matter with the expansion of the universe, we have that, 

until tea-, 



Pa 

p 



a a, 



(146) 



where the scale factor of the universe a oc 1/T. In view of 

the above, the requirement that axions do not overdose 
the universe translates into requiring that 



Pa 
P 



< 1, 



(147) 



cq 



since, by definition, peq ~ Pm, where pm is the energy 
density of matter. From Eqs. (145), (146) and (147), we 
obtain the bound 



< 



10^2 GeVy 



(148) 



where wc used that Teq « 2.8 x 10~^ GeV. Hence, an 
initial misalignment angle of order unity is possible only 
if Vo < 10^^ GeV. In view of Eq. (Ill), this, in turn, is 
possible only if n < 1. 

Thus, if n > 1 and wc insist on ~ Ij the only way 
that axion overproduction can be avoided is by consider- 
ing that the curvaton does dominate the universe before 
decaying. In this case, reheating due to curvaton decay 
dilutes the axion energy density because of a dramatic 
production of entropy. The entropy ratio at curvaton 
decay is easily estimated as 



«-> after 
re 



Pa_ 

Pi / REH 

9 0^vl 



"13/2 



Mi 



Mp 
W3/2 



(149) 
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where the subscript "REH" denotes the time of the cur- 
vaton decay, is the energy density of the background 
radiation due to the decay of the inflaton, and we have 
used Eqs. (Ill), (122), (128) and (135) considering also 
that fTosc ^ Ovq. The exact calculation multiplies [46] the 

result by a factor of 1.83 ^ 1, where 5* ^ 10 — 10^ is 

the effective number of relativistic degrees of freedom. As 
expected, the entropy production depends on the value of 
the coupling constant g. Hence, avoiding axion overpro- 
duction, which could overdose the universe, is expected 
to set another lower bound on g, more stringent than the 
one in Eq. (136) [47]. 

The fact that the case of curvaton domination requires 
a larger value of g [cf. Eq. (136)] is to be expected be- 
cause larger g implies that the inflaton decays earlier and, 
therefore, the energy density fraction Pa/p grows sub- 
stantially, allowing the curvaton to dominate the universe 
before its decay. The higher g is the more dominant the 
curvaton will be at its decay and, hence, the more diluted 
the axion energy density will become after the decay of 
the curvaton. 

A crucial further requirement for the dilution of the 

axion energy density is [30, 48 50] that the entropy re- 
lease occurs after the onset of axion oscillations. Hence, 
the requirement is 



> r, 



REH- 



(150) 



E. The evolution of the order parameter 

Let us now concentrate on the evolution of the order 
parameter v, which has to be such as to achieve the re- 
quired value for e. The order parameter is determined 
by the rolling \S\. When the cosmological scales exit the 
horizon, the field \S\ has to be slowly rolling because we 
need the order parameter to vary slowly enough not to 
destabilize the approximate scale invariance of the per- 
turbation spectrum [cf. Sec. IIIC]. Therefore, the Klein- 
Gordon equation for \S\ takes the form 



3if]5] mils'] ~0. 



(151) 



Using Eq. (104), the rate of growth of the order param- 
eter in this case can be easily found to be 



\S\ i""' \csH^ 



l\H. 



(152) 



The amplification factor can be found as follows. 
Using Eq. (9), we can write \S\ as a function of the num- 
ber N of the remaining e-foldings of inflation. Starting 
from Eq. (151) and after a little algebra, we obtain 



^ d\n\S\ e-2^«^- 



-2FsN 



cs dN 



1-e 



-2F,N 



(153) 



where corresponds to the phase transition which 
changes the sign of m|. Here, we used the fact that. 



by definition, 

\ml\^csHi:.csHl{l-e-'^^''^), 



(154) 



where = H{N^) and H^n = \A^/\/3Mp with Vm be- 
ing the scale of the inflaton potential as given in Eq. (7). 
Integrating Eq. (153), we get 



Imfm = (l-e--^-^)F-lnr^!!:^ 



(155) 



where [S*], = l^KA^,) and \S\^ = \S\{N^). 

Now, it is straightforward to check that, if 2FsN^. ^ 1, 
the RHS of the above equation tends to zero, which 
yields ]5]* « 15]x. This is also understood by observ- 
ing that, in this case, H'^{N) ~ Hl{l - e'^^^^) « 
Hl^, which means that S{H^)/H^ « 2FJNe-'^^'^ < 1, 
where S{H^) = - and 5N ^ N^- N^> 0. Thus, 
in the e- folding interval 5N , the effective mass-squared of 
S hardly changes: nigiN^,) « m|(A'x) = 0, i.e. the mass 
is very close to zero. This implies that \S\ remains frozen 
and, thus, \S\^ « \S\x- Moreover, the displacement of \S\ 
from the origin at the phase transition is determined by 
its quantum fluctuations, which means that 



27r' 



(156) 



For 2F,N^ » 1, iJ^ w (w H^) and, thus, \S\^ ~ 
if*/27r. So, under these circumstances, we have, for 
the amplification factor, e = £rnin) which is defined in 
Sec. niA. 
It is easily seen that, generally, 



£ = 



15)* ^ 15]* if* 
\S\o vo 



-H,. (157) 



Note that, since the required e is always much bigger than 

Emin, as wc saw from Eqs. (118) and (119), \S\^, 3> ff*. 
Furthermore, it is obvious from the above discussion that, 
for £ > Ernin, we need to have 



1 



2Nx 

Finally, when 2FsNy, <C 1, Eq. (155) reduces to 



cs \\SW 



+ iV*. 



(158) 



(159) 



In view of Eq. (152), the requirement in Eq. (43) takes 
the form 



^ {H^^ - Hi) « Hi 



(160) 



where we took into account Eq. (154). This inequality 
can be recast as 



3 I 1 - ' 



(161) 
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When 2FsN^ ^ 1, the above reduces to 
cs - 



< 1. 



(162) 



This means that, in this case, the cosmological scales 
must exit the horizon not much later than the phase tran- 
sition which changes the sign of ffig. 

Another issue to be addressed concerns the require- 
ment that 15*1 does slow roll at the time when the cos- 
mological scales exit the horizon, in contrast to the case 
when the slope of the potential is so small that the mo- 
tion of |5| is dominated by its quantum fluctuations. In- 
deed, since |to|(7J,)| <C H^, \S\ is effectively massless 
and, hence, it obtains a superhorizon spectrum of per- 
turbations of order H^,/2tt, much like a. In order for its 
quantum fluctuations not to dominate its motion, \S\ has 
to be outside the quantum diffusion zone. The condition 
for this to occur is 



dV 



d\S\ 



2\fnl{HMSU>H'. 



(163) 



Using Eqs. (104) and (154) and working as before, the 
above constraint is recast as 



In 



1^ 



> 2F,N^ - In ( — 



In 



1 



2 Vl 



o-2F,(iVx-A'.) 



(164) 



where we have also used Eq. (156). If 2FsN^ <C 1, the 
above equation becomes 



In 



In 



In 



7V| 



> 0. (165) 
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FIG. 2: Plot of the scalar potential VdS"! , |P|) which is de- 
fined in Eqs. (61)-(65) in units of M/^/Mp" with respect to 
15*1 and which are measured in units of Mi. Here, we 
have taken n = 5 and Mi = (mMp)^^^ . We have also cho- 
sen mp = niQ = — m| = m^, A = —9 m, Ao = Ai = A2 = 1, 
\P\ = IQI and 9s ~ 9p = 6q = so that the potential is min- 
imized according to Eqs. (87) and (95). The shifted val- 
ley is clearly visible. It starts at the value \P\ ~ 0.814 Mi 
when l^l = and arcs down towards the \S\ axis, pass- 
ing through the non-trivial minimum at |5|() ~ 0.683 Mi and 
|P|o ~ 0.772 Mi. The trivial valley at \P\ = can also be dis- 
cerned. 



The resulting scalar potential is shown in Fig. 2. 

Using Eq. (118), we obtain the value of the amplifica- 
tion factor necessary for the model to work: 



lO-^-^ridc 



(168) 



Now, let us assume, at first, that fidoc < !> he. the cur- 
vaton decays before domination. Then, Eq. (127) gives 



VIII. A CONCRETE EXAMPLE: n = 5 & 6> ~ 1 

From the bounds on n in Eqs. (124), (133) and (138) 
and also in view of Eq. (59), we see that not many choices 
for n are allowed. In fact, we can only accept the cases 
corresponding to I = 1,2 (i.e. n = 5,9) with the latter 
choice being marginal as far as the BBN constraint in 
Eq. (124) is concerned. Hence, to illustrate the above, we 
present an example taking n = 5 (i.e. / = 1 in Eq. (59)) 
and considering that the orthogonal axion assumes a ran- 
dom value after the phase transition, i.e. ~ 1. 

The bound in Eq. (124) suggests that this case is ac- 
ceptable provided that 



> 220 GeV. 



(166) 



The superpotential for the SM singlet superfields is com- 
prised only of the terms [cf. Eq. (60)] 



singlet 



[Ao^^ + \iS\PQf + \2{PQf]lMl. (167) 



10 



-4.75 



which means that 



5- 10"^-^ rid, 



(169) 



(170) 



This value lies well within the range given in Eq. (129). 

With these values, it is straightforward to show that 
all the relevant constraints are satisfied. For example, 
the constraint in Eq. (33), in the case when fidec < 1; 
becomes 



g fJdcc ^ {eC,f 



Mp 

"T-3/2 



n + 2 
n+1 



(171) 



which, for n = 5 and the value of g found above, can be 
easily checked to hold true. According to Eq. (130), the 
universe reheats at temperature 



10^ ridoc GeV, 



(172) 
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which is far from challenging the gravitino bound. 

However, the PQ symmetry breaking scale is found 
from Eq. (Ill) to be 

vo ~ lO^^-S GeV. (173) 

Such a high PQ scale results in axion overproduction un- 
less the original axion misalignment angle 0a is much less 
than unity. Indeed, in this case, Eq. (148) suggests that 

< 10"^ (174) 

should be enough to avoid axion overproduction. The 
above constraint is much less stringent than the con- 
straint on the present value of 6a, from CP violation in 
strong interactions (coming from experimental bounds on 
the electric dipole moment of the neutron): 0^ < 10~^, 
but it still undermines the motivation for the PQ symme- 
try, which is meant to explain how wc can get the present 
value of so small starting with a natural value of the 
initial 0a ~ 1- This is why it is preferable to consider the 
case when the curvaton decays after domination, where 
the axion energy density can be efSciently diluted by the 
entropy production from the curvaton decay. 

When the curvaton decays after domination, we have 
fidec ~ 1- Then, Eq. (168) becomes 

e~10-*-^ (175) 

Also, from Eq. (136), we obtain the bound 

g > 10"^-^ (176) 

As a result of the above, the reheat temperature after the 
end of inflation is found, from Eq. (130), to be 

Treh ~ g^ms/2Mp > 10^ GeV. (177) 

Satisfying the gravitino bound (Treh < 10^ GeV) sets a 
weak upper bound on g: g < 0.03. However, this bound 
is relaxed by the dilution of the gravitinos due to the 
entropy production by the curvaton decay. Indeed, in 
this case, according to Eq. (139), the hot big bang begins 
at the temperature 

Treh - 10 MeV, (178) 

which is close to the BBN bound, but does not violate it. 

Now, from Eq. (149), we find the entropy production 
ratio to be 

1^ ~ lO'^-^g. (179) 

The relic abundance of axions can be calculated by ap- 
plying the formulae of Ref. [43], where we take the QCD 
scale Aqcd = 200 MeV and ignore the uncertainties for 
simplicity. Comparing this relic abundance with the best- 
fit value of the cold dark matter (CDM) abundance in the 
universe from the measurements of WMAP [51], we find 
that, for Vo as large as the one shown in Eq. (173), the 



entropy of the universe must be increased by a factor of 
order 10^ at least. In the present model, this is actu- 
ally an overestimate since, as it turns out, the axions are 
generated, in the viable cases, after the domination of 
the oscillating curvaton field. Axion production in a uni- 
verse dominated by a coherently oscillating scalar field 
has been studied in Ref. [49]. In such a universe, ax- 
ion oscillations begin at a smaller temperature than in a 
radiation dominated universe. Consequently, the initial 
ratio of the axion number density to the photon number 
density is smaller. This then implies that, in this case, 
less entropy production is required. Taking this effect 
carefully into account, we obtain the bound 

g > 10-4 ■^ (180) 

which is more stringent than the bound in Eq. (176), 
as expected. A large coupling between the inflaton field 
and its decay products can be realized if the VEV of the 
s modulus corresponds to a point of enhanced symme- 
try. Note that, in this case, there is no moduli problem 
becaiisc! .s decays miich earlier than BBN. 

Moreover, we have to make sure that the entropy pro- 
duction occurs after the onset of the axion oscillations. 
The latter, according to Eq. (141), occurs at energy den- 
sity given by 

1 GeV, (181) 

where we used Eq. (173). Comparing the above equation 
with Eq. (178), we see that the requirement in Eq. (150) 
is satisfied. Thus, we conclude that axion overproduction 
can be avoided, due to entropy release at curvaton decay. 

It should be noted, however, that this mechanism of 
diluting the axions by the entropy produced when the 
curvaton decays after dominating the universe may lead 
[30] to a cosmological disaster. Generally, a sizable frac- 
tion of the curvaton's decay products consists of sparti- 
cles, which eventually turn into stable lightest sparticles 
(LSPs) in models (such as ours) with an unbroken mat- 
ter parity symmetry. The freeze-out temperature of the 
LSPs is typically much higher than the reheat temper- 
ature in Eq. (178). Actually, it is even higher than the 
energy density scale corresponding to the onset of ax- 
ion oscillations (see Eq. (181)). Consequently, the LSPs 
freeze out immediately after their production and can, 
subsequently, overdose the universe leading to a cosmo- 
logical catastrophe. 

In the present model, the PNGB curvaton can decay 
into a pair of squarks, sleptons, charginos or neutralinos 
with a decay width which can be comparable to the width 
of its main decay channel to a pair of Higgs particles (see 
Eq. (122)). We should, however, observe that our cur- 
vaton, being a SM singlet, couples to the charginos and 
neutralinos only via their Higgsino component. Taking fi 
to be much greater than the soft masses of the bino and 
wino, which is often encountered in various realizations 
of MSSM, we can ensure that the lighter charginos and 
neutralinos are predominantly gauginos. We can further 
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choose all the squark, slcpton and heavier chargino and 
ncutralino masses to exceed half the curvaton mass (see 
Eq. (166)) so that all the curvaton decay channels in- 
volving these particles are kinematically blocked. Thus, 
the only decay channels allowed are to a pair of lighter 
charginos or neutralinos. The corresponding rates can be 
easily suppressed by reducing the Higgsino component of 
these sparticles. Indeed, assuming that their Higgsino 
components are about 1%, we obtain a suppression fac- 
tor of order 10~*. In view of the fact that these sparti- 
cles have comparable masses to the Higgs particles, we 
then conclude that only a fraction of about 10~^ of the 
energy density of the universe soon after the curvaton 
decay consists of LSPs. This fraction is enhanced by a 
factor of about 10^ until the time teq of equal matter 
and radiation energy densities when the cosmic temper- 
ature Toq ~ 10~^ GeV and remains essentially constant 
thereafter. So, the cosmological disaster from the possi- 
ble overproduction of LSPs at the curvaton decay can be 
avoided. Moreover, the LSPs can contribute to the CDM 
in the universe. 

Using Eq. (119) with n = 5, we find 



10 



-12.5 



(182) 



Therefore, Eqs. (157) and (175) suggest that 

\S\^ - 10*H^. (183) 

The above can, in principle, be used in Eqs. (155) 
and (164) to constrain the parameters of the underlying 
model (e.g. Fg). 

A useful quantity to calculate in order to evaluate 
Eqs. (155) and (164) is the number of e-foldings which 
corresponds to the cosmological scales A''*. The cosmo- 
logical scales range from a few times the size of the hori- 
zon today <^ Hq^ down to scales ~ 10~^Hq^ correspond- 
ing to masses of order IQ^Mq with Mq being the solar 
mass. Typically, this spans about 13 e-foldings of in- 
flation. For the estimate of A^*, wc will choose a scale 
roughly in the middle of this range. More precisely, we 
will take the scale that re-enters the horizon at the time 
when structure formation begins, i.e. at the time t^q of 
equal matter and radiation energy densities. Then it is 
straightforward to obtain 



exp(7V*) 

T T 
JJ* ten 



inf 



vo 



(184) 



where we have used Eqs. (122) and (135) with CTosc ~ vq. 
Note that, remarkably, the above is independent from g. 
Putting if* ~ nia ^ 1713/2, we obtain 



38, 



(185) 



where wc have also taken Eq. (173) into accoimt. The 
number of e-foldings that corresponds to the decoupling 
of matter and radiation (when the CMBR is emitted) is 
roughly N^, + 1.5, while the one which corresponds to the 
present horizon is about iV« -|- 9. 



Using the above, let us attempt to investigate first the 
case when 



2FsN^ < 1. 



Substituting Eq. (159) in Eq. (162), we obtain 



(22.1 -mf) (f^-i) 



76 



1 



«1, 



(186) 



(187) 



where we also used Eqs. (156), (183) and (185), and the 

fact that ~ H^{1 — e~^^^^), which, in this case, is 

« 2FsNH?;,. In order to find how small the left hand 
^ 111 

side (LHS) of this inequality should actually be, we must 
observe that the contribution to the spectral index of 
density perturbations Ug originating from the evolution 
of V during inflation is —2H~^{v/v)^,, which is negative. 
Moreover, one can easily check that, in the present exam- 
ple, all the other contributions [52] to the spectral tilt for 
the curvaton are negligible. This is due to the fact that 
e 1 (see Eq. (175)) and, as it turns out, also c <C 1 
(see below). Using Eq. (152), one can further show that, 
under these circumstances, dris/dlnk is very small and, 
thus, no running of the spectral index is predicted in our 
model. For fixed Ug, the recent results of WMAP im- 
ply [51] that Us = 0.96 ± 0.02. Therefore, at 95% c.L, 
rig > 0.92. It is then obvious that the LHS of the in- 
equality in Eq. (187) should not exceed about 0.04. This 
requirement is met provided that 



N, 



r > 500, 



which, in view of Eq. (185), implies that 
iVx > 1.9 X 10^. 

Prom Eq. (159), we then obtain that 
cs < 2.4 X 10"''. 



(188) 



(189) 



(190) 



It can be checked that, with these values, the requirement 
in Eq. (165) is well satisfied. From Eqs. (186) and (189), 
one obtains 



Fs «: 2.63 X 10" 



c «: 7.89 X 10" 



(191) 



where we also used Eq. (8). Such a small c implies that 
modular inflation is not really of the fast-roll type and 
may last for a large number of e-foldings. Indeed, ac- 
cording to Eq. (13), iVtot > 10'^. In view of Eqs. (8) and 
(154), Eqs. (190) and (191) suggest that 

\ms\ < 1.55 X IQ-'^H^ and < 8.88 x lO'^H^ 

(192) 

with < 4.47 x 10 ^Hm- The above values for c and 
Cs are plausible (requiring only mildly tuned masses) 
but not very pleasing, since we would prefer c ~ cg ~ 1. 
(Note that larger values of result in more severe tun- 
ing of c and eg.) The small values obtained may be due to 
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the condition in Eq. (186), which we imposed to simphiy 
the problem. 

Therefore, let us consider, now, that 



2F,N^ < 1 and 2F«7Vx > 1. 
In this limit, Eq. (161) takes the form 
22.1 + ln(2^^,iV,) 



761n(2F,Af*) 



< 1. 



(193) 



(194) 



We find that the LHS of this inequality remains smaller 
than 0.04 provided that 2FsN^ < 0.004, which yields 



F, < 5.26 X 10"^ 



c < 1.58 X 10"^. 



(195) 



This is less fine-tuned than the values in Eq. (191). In 

order to estimate the lower bound on A'x corresponding 
to the upper limit on Fg, we approximate Eq. (161) for 
2FsN^ < 1, but any value of 2FsN^: 



[22.1 -ln(^ 



-)](l-e- 



76[(1 - e- 



ln(i 



2F^N. 



^) - 2F,iVJ 



< 1. (196) 



For 2FsiV* ~ 0.004, the LHS of this inequality is kept 
smaller than 0.04 even if 2FsA'x becomes as low as 4, 
which yields 

TVx > 3.8 X 10*. (197) 

Saturating the bounds in Eqs. (195) and (197), i.e. taking 

i?; ~ 5.26x10"^ (c~ 1.58x10"*) and N^:^ 3.8x10^, 

(198) 

we find from Eq. (155) that 



cs ^ 4.92 X 10" 



(199) 



which is a little more natural than the values in Eq. (190). 
It can be checked that, with the values in Eq. (198), the 
requirement in Eq. (164) is well satisfied. Using Eq. (13), 
it is easy to see that, in this case, the total number of 
e-foldings of inflation is 

iVtot ^ 6.6 x 10^ (200) 

In view of Eqs. (8) and (154), Eq. (198) suggests that 

\ms\^ 2.22 X 10-^ and ~ 1.26 x lO'^i?* 

(201) 

with ~ 6.32 x 10 ^ffmj which are not very different 
from the upper bounds in Eq. (192). After some inves- 
tigation, it can be realized that not much improvement 
can be made on the results. 

In conclusion, we see that, in the n = 5 case with 6* ~ 1, 
our model can work, typically, for values 



cs, c < 0(10-*) 
or, equivalently, for masses 

\ms\, rus < 0(10-2) if,. 



(202) 

(203) 



where ^ '713/2 ^ 1 TeV. Such values are quite natural 
and imply only a mild tuning on the masses of the rolling 
field 1 5*1 and the inflaton modulus. This is necessary be- 
cause the effective mass |tos| should remain small enough 
during the relevant part of inflation for l^l to be slowly 
rolling and the constraint in Eq. (43) to be satisfied. The 
condition for this is that cs or, equivalently, \ms \ be small 
[53]. Yet, a substantial variation of |5| from the phase 
transition until the time when the cosmological scales exit 
the horizon is necessary for obtaining the required value 
of the amplification factor (recall that the required e 
is always much bigger than Emin)- This is achieved with 
a small mass for the inflaton modulus, which leads to a 
large number of e-foldings. So, modular inflation cannot 
be of the fast-roll type in this case. The above flndings 
are similar to the ones in Ref. [18], despite the fact that 
there the example studied considered the case when the 
curvaton decayed before domination with n = 2. This 
suggests that the above results are quite robust. 

One may wonder why, since both the inflaton s and the 
held \S\ turn out to be light when the cosmological scales 
exit the inflationary horizon, we cannot use those flelds 
to generate the observed curvature perturbation. The 
reason is that, in contrast to the PNGB curvaton, the 
perturbations of those fields are not amplified. Hence, 
their contribution to the overall curvature perturbation 
is insignificant. Indeed, for the infiaton, we have 



1 



v:- 



5^377 \V'UM^ 57r 



^ ^^^jv.^ (204) 



Mp 



where the prime denotes derivative with respect to the in- 
flaton s and we have used Eqs. (6) and (8). For the above 
discussed values in Eqs. (185), (198) and (201) and for 
-ff* ~ ~ 1 TeV, Eq. (204) gives C ~ 10'^"^, which 

is much smaller than the observed value ^ ~ 2 x 10-^. 
Similarly, for |5|, we have 



^1^1 3 \s\ 



H* . 
eCr, 



(205) 



where \S\ = ||5| — \S\q\, 6\S\ the perturbation in \S\ and 
we have used the fact that |S'|o ^ t^o ^ jS'l* with \S\t, 
given by Eq. (183). For the values discussed above, 
C|s| ^ IQ-^^ <C C> where we considered that C ~ Co-- 



IX. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied modular inflation, which 
uses a string axion as the inflaton fleld. The inflationary 
scale (~ 10^°-^ GeV) is determined by the scale of gravity 
mediated soft SUSY breaking. Such low-scale inflation, 
even though it can still solve the flatness and horizon 
problems of the standard hot big bang cosmology, can- 
not generate the observed curvature perturbation (nec- 
essary to explain the CMBR anisotropy and structure 
formation) from the quantum fluctuations of the inflaton 
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field. However, we have shown that this type of modu- 
lar inflation can generate the appropriate amplitude of 
superhorizon curvature perturbation to account for the 
observations through the use of a suitable curvaton field. 

The curvaton field that we have used is a PNGB which 
is an angular degree of freedom orthogonal to the QCD 
axion field (that we called the orthogonal axion) in a class 
of SUSY PQ models. We considered models that gener- 
ate the PQ scale dynamically (by using flaton fields), 
while they also solve the problem of MSSM. In these 
models, one needs more than one SM singlet superfields 
to break the global U(1)pq symmetry (with the excep- 
tion of its matter parity subgroup). Hence, apart from 
the axion, there is at least one other angular degree of 
freedom (the orthogonal axion) which may be kept appro- 
priately light during inflation and can be responsible for 
the curvature perturbation in the universe. This could be 
achieved if the potential possesses a valley of minima with 
a negative inclination and the system happens to slowly 
roll down this valley during the relevant part of inflation 
with some of the SM singlet fields acquiring values much 
smaller than their vacuum values. Under these circum- 
stances, the orthogonal axion may be kept light during 
inflation and its perturbation from inflation may be later 
amplifled as the SM singlets acquire their vacuum val- 
ues accounting for the observed curvature perturbation 
in the universe. 

Following this promising idea, we have attempted to 
construct appropriate curvaton models using two SM sin- 
glet superfields P and Q, charged under the PQ symme- 
try. However, we have shown that it is not possible to 
construct suitable curvaton models by using only two SM 
singlet superfields, because, in this case, the orthogonal 
axion mode cannot avoid being massive during the rel- 
evant part of inflation. Hence, we studied PQ models 
which involve a third superfleld S and specifled the gen- 
eral conditions under which these models can contain a 
suitable PNGB curvaton. Actually, we have shown that 
they possess a shifted valley of minima at almost con- 
stant values of |P| and \Q\ of the order of their vacuum 
values. The value of |S'|, however, which parametrizes the 
valley, can be kept much smaller than its vacuum value 
during the relevant part of inflation. Also, there exists 
an orthogonal axion which remains light during inflation 
and can serve as curvaton. 

For definiteness, we considered a concrete class of mod- 
els of this category where the superfield S has vanishing 
PQ charge (see also Ref. [10]). These models are simple 
extensions of MSSM based on the SM gauge group and 
possessing, besides the global anomalous PQ symmetry, 
a global U(l) R symmetry and a discrete Z2 symme- 
try. The superpotential of the models includes all the 
terms satisfying these symmetries. The baryon and lep- 
ton numbers are automatically conserved to all orders in 
perturbation theory as a consequence of the R (and the 
PQ) symmetry. 

To stiidy the cosmology, we have also taken into ac- 
count that the SM singlet fields in these models {P, Q 



and S) are expected to receive SUGRA corrections to 
their soft SUSY-breaking masses-squared (as well as their 
soft SUSY-breaking trilinear A-terms) of order set by 
the Hubble parameter. We consider the corrections to 
the masses-squared to be positive in all cases. Since our 
infiation model has H^, of order the electroweak scale, 
which is also the scale of the soft masses, these SUGRA 
corrections do not seriously affect the physics with the ex- 
ception of the S field, whose soft mass-squared is taken 
to be negative. This is intentional in order to facilitate a 
phase transition during inflation which sends \S\ rolling 
away from the origin along the shifted valley and reaching 
its true minimum by the end of inflation (cf. Ref. [18]). 
As a result, the order parameter of our PNGB curvaton 
field increases substantially after the cosmological scales 
exit the inflationary horizon. This amplifles the curva- 
ton perturbation according to the mechanism presented 
in Ref. [13] and the observed value of the curvature per- 
turbation in the universe can be achieved despite the low 
inflation energy scale. 

We have investigated in detail the above scenario and 
showed that the requirements for a successful curvaton 
put important constraints both on the choice of model 
and also on the model parameters. Indeed, we have 
shown that only a few members of our class of PQ models 
are eligible for successful curvaton. We then have con- 
centrated on a particular such model and, using natural 
values for the model parameters, we have studied analyt- 
ically its performance as curvaton modc^l. We found that 
the model can indeed work successhilly in the context of 
modular inflation with only a mild tuning of the inflaton's 
mass and the mass of S: \ms\,ms < O.OliJ*. The bound 
on jmsj comes from the requirement that |5| be slowly 
rolling when the cosmological scales exit the horizon, oth- 
erwise the scale invariance of the spectrum of curvature 
perturbations will be destabilized. The bound on rus, 
on the other hand, originates from the large variation 
of \S\ between the phase transition and the exit of the 
cosmological scales from the horizon, which is necessary 
for obtaining the correct amplification of the curvaton 
perturbation. This bound implies that modular infiation 
cannot be of the fast-roll type and may last for a large 
number of c-foldings. 

In this model, the PQ scale turns out to be quite large 
(it is comparable to the scale of grand unification). It 
is actually well above the standard cosmological bound 
from the requirement that the universe is not overdosed 
by axion overproduction. However, overclosure of the 
universe can be avoided, in this model, by diluting the 
primordial axions through the entropy release by the de- 
cay of the curvaton field. For this to be effective, the 
infiaton modulus has to decay early enough so that the 
curvaton may well dominate the radiation background. 
Consequently, we need a comparatively large decay cou- 
pling constant for the infiaton, which is possible if the 
VEV of the infiaton corresponds to an enhanced symme- 
try point. 

After the curvaton decays, the hot big bang begins. 
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This occurs not long before BBN. So, baryogenesis has 
to take place soon after the time of the curvaton decay 
at the latest. Moreover, it is clear that, in the present 
case where the curvaton decays after dominating the en- 
ergy density of the universe, baryons must be generated 
through the decay of the curvaton (or of some of its decay 
products) , since otherwise there will be [25] an unaccept- 
ably large baryon isocurvature perturbation. It is also 
obvious that leptogenesis [54] (or electroweak baryogen- 
esis) cannot work here since the reheat temperature is 
very low for the non-perturbative electroweak sphaleron 
effects to operate. Therefore, almost the only viable op- 
tion is that the observed baryon asymmetry of the uni- 
verse is directly generated by the decay of the curvaton 
(or of its decay products). One possibility is that the cur- 
vaton has suitable baryon (and lepton) number violating 
decay channels via non-renormalizable Lagrangian oper- 
ators of higher order with decay widths comparable to 
its main decay width to Higgs particles. This, of course, 
requires an appropriate extension of our model, which. 



as it stands, has exact baryon (and lepton) number con- 
servation. It may be possible [55] to obtain the required 
Lagrangian operators by embedding our model in a larger 
scheme with (large) extra dimensions. This baryogenesis 
issue deserves further study, which we postpone for the 
future. 

The curvaton model analyzed in this paper can be con- 
sidered to accommodate [56] low-scale inflationary mod- 
els other than modular inflation of the type considered 
here. However, modular inflation due to string axions is 
one of the better theoretically motivated low-scale infla- 
tion models in the literature. 
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